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o : 

Abstract. In this article we consider the af finely-rigid body moving in the three-dimen- 
sional physical space and subject to the Kirchhoff— Love constraints, i.e., while it deforms 
homogeneously in the two-dimensional central plane of the body it simultaneously performs 
one-dimensional oscillations orthogonal to this central plane. For the polar decomposition 
OO ' we obtain the stationary ellipsoids as special solutions of the general, strongly nonlinear 

equations of motion. It is also shown that these solutions are conceptually different from 
pH ' those obtained earlier for the two-polar (singular value) decomposition. 

Key words: affinely-rigid bodies with degenerate dimension; Kirchhoff-Love constraints; 
_C ■ polar decomposition; Green deformation tensor; deformation invariants; stationary ellipsoids 
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oo ' 1 Introduction 

-xf : 

The special interest in the present work is devoted to the classical description of an affinely-rigid 
(homogeneously deforming) mechanical system subject to the Kirchhoff-Love constraints. We 
know that the standard continuum theory as well as some fundamental theories deal with such 
\ objects as membranes, plates, discs, etc. So, the main contribution of this work is to present 

a toy model for the analytical description of the above-mentioned objects. 

The structure of this article is as follows: Firstly, we will present the main notions about the 
concept of the affinely-rigid body, as a generalization of the metrically-rigid one, and of its special 
case, i.e., the affinely-rigid body with degenerate dimension. Secondly, for convenience of the 
Reader the main results obtained earlier for the case of two-polar (singular value) decomposition 
are remembered. Thirdly, an alternative (polar) decomposition is introduced and the equations 
of motion for our toy model are obtained for the general form of the inertial tensor, i.e., when 
^1/^2/ ^3' And finally, three main branches of special solutions (stationary ellipsoids) for our 
strongly nonlinear equations of motion are gathered in the form of Proposition [TJ Additionally 
some remarks about the complementarity of the obtained results to those described in our 
previous work [6] are presented in the Summary. 

So, let us remind some basic facts generally concerning the notion of affinely-rigid bodies [141 

EH EH EE]. 

Let (M, V, — >) be an affine space and (M, V, — >, g) be the corresponding Euclidean one, where 
M is a physical space in which the classical system of material points (discrete or continuous) 
is placed, V is a linear space of translations (free vectors) in M, and g € V* <g> V* is the metric 
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tensor. Also let us introduce an affine (N, U, — >) and the corresponding Euclidean (N, U, — rf) 
spaces, where N is the material space of labels which are assigned to every material point of our 
body in some way, U is the corresponding linear space of translations in N, and rj £ U* <8> U* 
is the metric tensor. Then the position of the o-th material point at the time instant t will be 
denoted by x(t, a) (x € M, a G N) and an affine mapping from the material space into the 
physical one is as follows: 

x i (t,a) = r i (t) + ip i A (t)a A , 

where <p(t) is a linear part of the affine mapping (ip is non-singular for any time instant t), i.e., 
(p(t) £ IA(U,V), where LI (£7, V) is a manifold of linear isomorphisms from the linear space U 
into the linear space V, r(t) is the radius- vector of the centre of mass of our body if in the 
material space the position of the centre of mass is a A = 0. If the system is continuous, then the 
label a becomes the Lagrangian radius-vector (material variables) and x becomes the Eulerian 
radius- vector (physical variables). Thus, at any fixed t 6 I the configuration space Q of our 
problem is given by the following expression: 

Q = Afl(/V, M) = Q tr x Q int = M x L1(U, V), 

where "tr" and "int" refer to the translational (spatial translations) and internal (rotations and 
homogeneous deformations) motions respectively. 

The considered system is called an af finely-rigid body [HI [T5], HU H7J Q2], i.e., during any 
admissible motion all affine relations between constituents of the body are invariant (the ma- 
terial straight lines remain straight lines, their parallelism is conserved, and all mutual ratios 
of segments placed on the same straight lines are constant). The conception of the af finely- 
rigid body is a generalization of the usual metrically-rigid body, in which during any admissible 
motion all distances (metric relations) between constituents of the body are constant (see, for 
example, [H H]). 

In this article we concentrate mainly on the case of such an affinely-rigid body that is subject 
to the additional constraints, i.e., it can deform homogeneously in the two-dimensional central 
plane of the body and simultaneously performs one-dimensional oscillations orthogonal to this 
central plane. Then the material space ./V is presented as the Cartesian product R + x R 2 and 
the group of material transformations has the form R + x GL (2,]R), where R + is the dilatation 
group in the third dimension and the material transformations in R 2 act as in the case of the 
usual affinely-rigid body with degenerate dimension 112], 

We can identify configurations <I> : R 3 — > R 3 with the pairs (g,^p), where ip describes the 
immersion of the central plane in the physical space, i.e., analytically ip 1 a is the 3x2 matrix. 
An element (k, B) acts on (g, ip) as follows: 

(k, B)eR + x GL(2, R) : (g, ip) h-> (kg, tpB). 

The conservation of orthogonality of the direction of dilatations to the central plane means that 
the matrix 



$ 2 1 $ 2 9 $ 2 



3 



1 * 2 * 3 

3 



<3? 3 i < 3? 3 2 $ 3 



fulfils the condition that third column has to be proportional to the vector product of first 
and second ones. If we consider 3> a i, & b 2, a,b = 1,2,3, as independent and arbitrary, then 
<£ a 3 = £ e a bc& b i& c 2, where e a bc is the completely antisymmetrical Levi-Civita (permutation) 
symbol, £ is the parameter which depends both on the variable describing one-dimensional 
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oscillations orthogonal to the central plane of the body and on the ones describing the state of 
deformation in this central plane, e.g., for the two-polar (singular value) (HD and the polar Q 
decompositions we have respectively that 

nwo- polar ' ^polar ^ q/2 ' 

where the meaning of variables A, /i, a, £, £, g is clear from the expressions © and © below. 

The above-described orthogonality is well known in the theory of plates and shells as the 
Kirchhoff-Love condition [8]. 



2 Two-polar decomposition 

In [B] we discussed the language of the two-polar (singular value) decomposition: 
$ (jfe; A, n, g;6) =R (3b) D (A, fi, g) U (6») _1 , A, fi, g > 0, 



1 



where R,U € SO(3,M) are proper orthogonal matrices (whereas A; is a rotation vector, i.e., 
a non-normalized vector codirectional with the rotation axis whose magnitude is equal to the 
rotation angle) and D is diagonal, i.e., 



D(\,fj,,g) 



A 
n 
g 



cos 9 sin 9 
— sin 9 cos 9 
1 



(2) 



Then the co-moving angular velocities for R- and [/-tops [17} [18] are as follows: 
oj = R^R = R T R = 



co 3 -uj 2 
—UJ3 uj\ 
0J2 —uii 



-UJ, 



and 



For 6 and <I> T we have the following expressions: 



0-10 
1 




i) 1 



-i). 



$ = R(D + ujD - DtyU' 1 , 



T 



U(D + $D - Du)R T . 



The kinetic energy is assumed to have the usual form (we have only to substitute the con- 
straints): 

T= ±Tr(j$ T <i») = \Ti{U~ l JU[D + -dD - Duj] [D + ujD - D-d]), 

where J € U <8> U is the twice contravariant, symmetric, non-singular, positively-definite tensor 
describing the inertial properties of our affinely-rigid body. If we take J in the diagonal form 
J = Diag (Ji, J2, J3), then the above kinetic energy can be rewritten as follows: 



T 



J\ cos 2 9 + J 2 sin 2 9 ( d\\ J\ sin 2 + J 2 cos 2 9 ( d\x 



dt + 



dt 



+ 



J3 fdA 2 (Ji sin 2 6> + J 2 cos 2 0) /i 2 + J 3 g 2 



2 \dt J 
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+ 



[Jx cos 2 9 + J 2 sin 2 9) A 2 + J 3 g 2 



w| + ( J\ + J 2 ) A/LiW; 



d<9 
eft 



+ (Ji - J 2 )sm29 



d[i dX \ d9 



dt 



dt J dt 



dfi dX\ 



A* "77 - A— — + A— - fl— W 3 + \fXUJiLU2 



dt 



dt J 



+ 



+ 



( Ji cos 2 + J 2 sin 2 0) A 2 + ( Ji sin 2 + J 2 cos 2 0) ^ 2 

2 

( Ji sin 2 9 + J 2 cos 2 0) A 2 + ( Ji cos 2 9 + J 2 sin 2 0) ^ 2 



d0 



(3) 



The above expressions significantly simplify when we consider the isotropic case in the central 
plane of the body, i.e., when we have J\ = J 2 = J. Then 



J 
2 



dX 
~d~t 



+ 



d/i 
~dt 



+ 



J 3 / dg 



dt 



+ 



Jfl 2 + J 3 Q 2 2 



+ ^^- 2 2 + 2JA^ 3 | + 7j(A 2 + A i 2 ) 



"3 + 



We also remind here that the corresponding expression for the kinetic energy in the canonical 
variables has the following form: 



r 



+ 



2 ( J/i 2 + J 3 f? 2 ) 2 ( ja 2 + J 3 Q 



(A 2 + ^ 2 ) (s 2 + p 2 ) - 4A/ip e g 3 pj+p 2 , ^ 



2 J (A 2 



2J 



2J 3 



Then introducing some modelled potentials in [6] we obtained the Hamiltonian (total energy) 
and calculated the corresponding equations of motion for the isotropic case with the help of the 
Poisson brackets. In the present article we concentrate mainly on the alternative decomposition, 
i.e., the polar one. The main advantages of this decomposition are the more physically intuitive 
division on three main terms in the kinetic energy expression (see the formulas below) 
and the possibility to obtain the equations of motion in the quite simple form (see the expressions 
(]lip - (jl7p below) even for the general case, when the inertial tensor is not isotropic in the central 
plane (Ji ^ J 2 ). 



3 Polar decomposition 

Instead of ([1]) we can also use the language of the polar decomposition, i.e., 
$ (7c; a, £, (,g) = L (7c) 5 (a, f , C, g) , 



(4) 



where L 6 SO(3, R) is a proper orthogonal matrix and S G Sym(3,R) is symmetrical. The 
connection between the polar and two-polar decompositions is given by the following expressions: 

L = RU-\ 








V 3 -V2 




= -u T = 


-u 3 


Vl 


= U(u-d) U^ 1 




V2 


-vi 









uj 3 + 9 


—oji sin 9 — 1^2 cos 9 


-U3 - 9 







u\ cos 9 — L02 sin 9 



coi sin 9 + L02 cos 9 w 2 sin ( 



UJ\ cos 1 
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,5' 



£ a 
a C 
Q 



UDU' 



A cos 2 9 + [i sin 2 9 (A — /x) sin 9 cos 6 
(A — fi) sin 9 cos 9 A sin 2 9 + /x cos 2 




(5) 



and then the Green deformation tensor, which is not sensitive with respect to the left orthogonal 
mappings, is as follows: 



G = = S 2 



e 2 + a 2 (£ + C) a 
(C + a ( 2 + a 2 











UD 2 U^ 



A 2 cos 2 (9 + ii 2 sin 2 (9 (A 2 - fi 2 ) sin (9 cos 9 
(A 2 -/i 2 ) sin 9 cos A 2 sin 2 9 + fi 2 cos 2 6> 
q 2 

where for the positive definiteness the parameters have to fulfil the conditions 

£ = A cos 2 9 + fi sin 2 9 > 0, C = A sin 2 + /i cos 2 6» > 0, 



£C - « 2 = A/U > 0, 



g > 0. 



For the polar decomposition we can as well introduce the concept of deformation invariants fC a , 
a = 1,2,3, which may be chosen, e.g., as the eigenvalues of the symmetric matrix G: 

det [G - KI 3 ) = 0, 

where I3 is the 3x3 identity matrix, and the solutions are as follows: 



/Ci, 2 = h(f + ( 2 + 2« 2 ± (£ + C) V(£ " O 2 + 4a 2 , 



The above deformation invariants are not sensitive with respect to both the spatial and material 
rigid rotations (isometries). 

Let us consider the Lagrangian L = T — V (<&) and then the Hamiltonian H = T + V ($>), 
where the kinetic energy ([3]) can be rewritten for the polar decomposition as follows: 



T — T rot + T rot _def + Tdef ' 



(6) 



where 



'rot 



Jia 2 + J 2 C 2 + J3Q 2 2 J1Z 2 + J 2 u 2 + J3Q 2 2 
v{ H Uo 



+ 



2 1 2 

Jit 2 + J2C 2 + (Ji + J2) a 2 2 



(Ji£ + J 2 C) az/if 2 



(7) 



describes the coupling between the angular velocity v of the L-top and deformation matrix S, 



-rot-dcf 



Jia Tt- j2a ^- {J ^- j2C) nr 3 



describes the connection between the angular and deformation velocities, and finally 



Ji + J2 ( da\ 2 J\ ( d£\ 2 J2 ( d(,\ 2 J3 / dg^ 



tdcf 



dt ) + 2 



dt 



+ 



r/f 



+ 



dt 



(8) 



(9) 



describes the kinetic energy of the deformation oscillations, whereas the potential term V ($) 
depends on $ only through the Green deformation tensor G = S 2 , i.e., the potential term 
adapted to the polar decomposition is a function only of a, £, £, and g. 
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Performing the Legendre transformation we obtain that 
= (Jia 2 + J 2 C 2 + hQ 2 ) v x - (Jig + J2C) «^2, 
= (Jig 2 + J20? + ho 2 ) v 2 - {Jxi + J2C) atvi, 
~- ( Jig 2 + J2C 2 + (Ji + J2) « 2 ) ^3 + Jiai - J 2 ag - (Jig ~ J2O 
: (Ji + J 2 )d-(Jig- J 2 g)i/ 3 , 
Ji (g + 0^3) , 

J 2 (g - a^J > 

= J3P, 

where 7Tj are canonical "spin" variables conjugate to angular velocities V{. 
Therefore after inverting the above dependencies, i.e., 

= (Jig 2 + J2O? + J3Q 2 ) 7Tl + (Jig + J 2 Q a7T 2 

Ul ~ J X J 2 {a 2 - gC) 2 + [Jig 2 + J2C 2 + (Jl + Ji) « 2 ] J3Q 2 + J 3 V 
(Jig + J 2 g) avr! + {.ha 2 + J 2 C 2 + J3Q 2 ) tv 2 





dT 


TV 1 = 


dv\ 




dT 


7T 2 = 


dv 2 




dT 


" 3 


dv 3 




Ol 


Pa = 


da 




dT 


= 


dg 




dT 


PC = 


~d~c 




dT 


Pe = 


dg 



a.. 



u 2 



J1J2 (a 2 - gg) 2 + [Jig 2 + J 2 C 2 + (Jl + Ji) a 2 ] J3Q 2 + J 3 V 
(Jl + J 2 ) fa + a (p ( -Pt)] + (Jig - J 2 g)Pa 



JiJ 2 (g + C) 2 

da = (Jig - J 2 Q fa + a (gc ~ gg)j + (Jig 2 + J2C 2 ) Pa 

dt JiJ 2 (g + C) 2 

dg _ P£ _ ( Jl + J 2 ) fa + a (P( - gg)] + (Jig - J 2 g) Pa 

rf* ~ Ji ° JiJ 2 (g + C) 2 

Jg _ PC _l_ q ( Ji + J 2 ) fa + a (P{ ~ Pt)] + (Jig ~ J2C) Pa 
dt~J 2 a j lJ2 ^ + Q 2 

dg_ _Pe_ 
dt~ J 3 

we obtain the kinetic energy in the canonical variables as follows: 
(Jig 2 + J 2 a 2 + J 3 g 2 ) tv 2 + (J ia 2 + J 2 g 2 + J 3 Q 2 ) TV 2 



T 



2 (Ji J 2 (a 2 - gg) 2 + [Jig 2 + J 2 g 2 + (Ji + J 2 ) « 2 ] J3Q 2 + J 3 V 

(Jig + J 2 Q aTV\TV2 

Ji J 2 (a 2 - gg) 2 + [Jig 2 + J 2 C 2 + (Ji + J2) a 2 ] J 3 g 2 + J 2 g* 



+ 



Jl + J 2 r , , o2 , Jlg 2 + J 2 C 2 2 

H ; 3 fa + « (Pc ~~ Pus H ; 

2JiJ 2 (g + g) 2 2J!J 2 (g + C) 2 

From the above kinetic energy expressions dS])-© one can see that the generalized velocities a, 
g, g corresponding to a, g, g and other variables describing the motion in the central plane of the 
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body are separated from the generalized velocity g describing the one-dimensional oscillations 
orthogonal to this central plane. The same can be said also about the above expression in 
the canonical variables (jlOp . i.e., the momentum p s conjugated to g is orthogonal (in the sense 
of metrics encoded in the kinetic energy expression) to the other canonical momenta. Hence, 
the most simple are those dynamical models in which also the isotropic potential will have the 
separated form: 

V (a, £, C, Q) = Vpi anc (a, £, C) + V e (g) , 

where as the potential V g we can take, e.g., the following potential which describes the nonlinear 
oscillations and is in accordance with the main demands of the elasticity theory, i.e., 

V 8 {q) = - + \q 2 , a,b>0, 
g z 

where the first term prevents from the unlimited compressing of the body, whereas the second 
one restricts the motion for large values of g, i.e., prevents from the non-physical unlimited 
stretching of the body. 

So, the Hamiltonian (total energy) can be written as follows: 

H = T +Vp Lam (a,Z,Q + V e (g), 

where T is taken in the form of (jlOp . Then the equations of motion can be calculated with the 
help of the following Poisson brackets: 

The only non-zero basic Poisson brackets are 

{a,p a } = = {C,Z><} = {QiPq) = !) = -eij k 7r k , 

where the former expressions follow directly from the definition of the Poisson bracket and the 
latter ones are based on the structure constants of the special orthogonal group SO (3,R). 
First of all, let us rewrite the kinetic energy (jlOp in a more symbolic way, i.e., 

n(7Ti,7r 2 ) T(7r 3 + a(p c -pg),p a ) Pg_ jg_ _pj 
2S 2J 1 J 2 {C + Cf 2J i 2J 2 2J s' 

where 

S = Ji J 2 (a 2 - CCf + [Ji? + «/ 2 C 2 + (Ji + J2) a 2 ] J 3 g 2 + jjg\ 
and two expressions built of the canonical momenta are as follows: 

Q (vri, tt 2 ) = ( JiC 2 + J 2 a 2 + hp 2 ) vr? + 2 ( Jtf + J 2 C) avrivr 2 

+ (Jia 2 + J 2 ( 2 + hp 2 ) ttI 
T (7T 3 + a (p c - p € ) ,p a ) = (Ji + J 2 ) [tt 3 + a (p c - p f )] 2 + {J^ 2 + J 2 C 2 ) p 2 

+ 2 (Ji^ - J 2 C) [vr 3 + a fa - p{)] p a . 

Then we obtain the following equations of motion: 

dTT! _ [(M + h() QTTi + ( JiQ 2 + J 2 C 2 + J 3 Q 2 ) 7T 2 ] 7T 3 

dt E 
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vr 2 [( Ji + J 2 ) [vr 3 + a (p c - )] + (Jig - J 2 C) Pa] 



dg s 

J 2 CP a _ J 2 [tT3 + a (PC _ Pl>] Pa . T (^3 + OL (p ( - pg) , p a ) 



ft (7Ti,7r 2 ) 



(11) 

(12) 
(13) 



dTT 2 _ [( Jlj 2 + J 2 Q 2 + hQ 2 ) 7T1 + (Jig + J 2 Q «7T2] ^3 

dt ~ S 

_ 7Ti [(J! + J 2 ) [7T 3 + a (p ( - pg)] + (J^ - J 2 C) p a ] 

JiJ 2 (g + C) 2 

^3 _ (Jig + J 2 Q g (A - *D + [ J l (" 2 - ^) + J 2 (C 2 - « 2 )] ^2 

dt S : 

da 3Vpl anc (J 2 7T^ + Ji7T^) a + (Jig + J 2 C) VTi7T 2 

(it 3a 3 

2JiJ 2 a(a 2 -gC) + (Ji + J 2 )aJ 3 £ 2 

H ^ ^ (VT1,7T 2 ) 

_ ( Ji + J2) [tt 3 + a (pc ~ Pg)] + (Jig - J 2 Q Pa / _ , (u , 

JiJ 2 (g + C) 2 c n 

dg 3y p i anc Jig^ 2 + Jiavrivra JiJ 2 C(gC-a 2 ) + JigJ3^ 2 , s 

Tt = h + s 2 " n I* 1 '*** 

_ JlgPq + Jl [^3 + « (PC - Pg)] Pa T(7T 3 + a(p C -pg),p a ) 

JiJ 2 (g + C) 2 JiJ 2 (g + C) 3 

dC _ deplane J2CA + J 2 Q7Ti7T 2 JlgJ 2 (gC ~ « 2 ) + JjgJsg 2 



-+ v °' v ^ (16) 

JiJ 2 (g + C) 2 JiJ 2 (g + C) 3 

^ - ^ (tt? + vr 2 ) + || [Jig 2 + J 2 C 2 + ( Ji + J 2 ) a 2 + 2 J 3P 2 ] fa,**) • (17) 



dt (if? 

The structure of the above expressions implies that even in the simplest case of the completely 
separated potential the dynamical coupling between the parameter describing one-dimensional 
oscillations orthogonal to the central plane of the body and the variables living in this central 
plane is present. 



4 Stationary ellipsoids as special solutions 

Our equations of motion (llip - (|17p are strongly nonlinear and in a general case there is hardly 
a hope to solve them analytically. Nevertheless, there exists a way for imaging some features 
of the phase portrait of such a dynamical system, i.e., we have to find some special solutions, 
namely, the stationary ellipsoids \14\ I15|. which are analogous to the ellipsoidal figures of equi- 
librium well known in astro- [2j and geophysics, e.g., in the theory of the Earth's shape [3]. 

In the case of the two-polar (singular value) decomposition (JTJ) we obtained the above- 
mentioned special solutions just putting the deformation invariants A, p, g and the angular 
velocities u, 1? equal to some constant values [6]. But now, in the case of the polar decomposi- 
tion (j3J), we see that the Green deformation tensor G, therefore the deformation matrix S, and 
the angular velocity v of the L-top have to be constant [IS], i.e., 

dt dt y ' dt K 1 dt dt y ' 

This means that the L-top performs the stationary rotation, i.e., if at the initial time t = we 
have that the configuration of the body is Lq, then at the time instant t the configuration will 
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be as follows: 
Looe", 

where o is the function composition symbol. We see that the whole affinely-rigid body, which 
at the initial time t = has the internal configuration $0 = ^0° S, at the time instant t will be 
in the following configuration: 

$(t) = L oe' /t oS = e 9t oL oS = e m o $ , (18) 

where v = Lq o v o Lq l . 

Proposition 1. While the affinely-rigid body rotates in the stationary way around the axis fixed 
in the physical and material spaces, the deformation and the angular velocity of rotation are not 
independent and related by some algebraic expressions. 

Proof. The trajectories of type (|18|) are the orbits of the Euler (spatial) action of the one- 
parameter orthogonal group \e vt : t 6 R} C SO (V, g). Nevertheless, during the motion the 
configurations of the body are deformed. At the same time the Green deformation tensor does 
not perform any oscillations, it is constant. This means that an equilibrium is set between the 
centrifugal forces coming from the rotation of the body and the elastic forces coming from the 
fact that S ^ Idf/. This kind of equilibrium is possible only if between the constant values 
of v, S are set some algebraic relations that guarantee the balance of the above-described forces. 
These algebraic relations between v, S obviously come from the equations of motion (|lip -(|17 p . 
Hence, we propose to divide them into the following three main branches: 

(i) v\ / 0, U2 = V3 = 0, then m, tt 2 i=- and ir 3 = p a = p^ = p^ = p e = 0; 

(ii) v 2 ^ 0, v\ = v 3 = 0, then 7Ti, 7r 2 7^ and = p a = P£ = PC = Pg = 0. 
For the first two cases the relations take the same form, i.e., 

dVplanc _ ( J 2^1 + JlA) a + (M + J2C) 

da S 

2 J x J 2 a (a 2 - fC) + (J1 + J2) a J3Q 2 „ , . /in , 
H %2 i2(7ri,7r 2 ), (19) 

S^plane Jlfrj + J X a^ 2 JlJ 2 ( (fC ~ « 2 ) + M^Q 2 „ , , 

— -Qp — = n 1 w2 S2(vri,7r 2 ), (20) 

dVplane + J2aiTl1T 2 J^J 2 (£C ~ Oi 2 ) + J 2 (J 3 2 

^ = 5 + =^ fi(7Tl,7r 2 ), (21) 

^ = ~ W + rt) + J ~§ [ J ^ 2 + J < 2 + ( J i + J 2) « 2 + 2 -hp 2 ] n K *a) , (22) 

with the compatibility condition 

(Jit + J2O a (tt? - tt|) + [Ji (a 2 - e 2 ) + J2 (C 2 - a 2 )] 7r l7 r 2 = 0. 

We see that, while our parameters v\ or V2 take completely arbitrary constant values, the above 
equations (|19p - ([22p describe their interrelation with the elements of the symmetrical matrix S, 
i.e., with q, £, £, g. 

(ra) 1/3 7^ 0, z^i = z/ 2 = 0, then TT3,p a ,P£,P( 7^ and tv\ = 1^2 = Pg = 0, whereas 

(J2C ~ JlO ^3 



Ji£, 2 + J2C 2 + {Ji + J2)a 2 ' 
JiC + J2C 2 + (Ji + J2) a 2 ' 
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J 2 a7r 3 

P< ~ ~ Jie + .h( 2 + (Ji + J2)a 2 ' 
So, for the third case we obtain the following relations: 
dVpianc _ ( Ji + J2) [vr 3 + a (p c - p f )] + ( Ji£ - J 2 C) 



da j lJ2 (£ + C )2 

deplane = JigPq + J\ [tt 3 + a (gg - g|)] p a T (tt 3 + a (p c - pg) , p a ) 

^ j!j 2 (e + o 2 JiJ 2 (e + o 3 

deplane _ ^CPq ~ ^2 [tT 3 + «(?>£- ^)] Pa . T (vr 3 + a -p^),P 



(Pt-Pd, (23) 



(24) 
(25) 



d( jij 2 (e + o 2 (f + o 

0. (26) 



This time our parameter f 3 has a completely arbitrary constant value and the above equations 
(|23l) - ([26]l describe the way in which a, £, £, Q are related to it. ■ 



Remark 1. It should be mentioned that the name "stationary ellipsoids" is not the most 
adequate for the description of the above-obtained stationary solutions. Of course, for the 
non-restricted affinely-rigid body we can visualize this kind of special solutions as follows: 

• At the beginning the body stays in the equilibrium configuration. 

• Then we switch on some mechanical device which deforms our body in the homogeneous 
way, i.e., this deformation is the superposition of three material stretchings with the 
coefficients D\, D 2 , _D 3 , where Di, i = 1,2,3, are the diagonal elements of the defor- 
mation matrix D in the two-polar decomposition $ = RDU^ 1 . 

• After this we start to rotate our mechanical device with the constant angular velocity i? 
around one of the main axes of the Green deformation tensor G = <E> T <I> = UD 2 U~ 1 so 
that the state of material deformation follows this movement of the device with the same 
angular velocity (note that our body itself does not rotate!). 

• And finally, the whole system consisting of the already rotating mechanical device and the 
body starts also to rotate with the constant angular velocity u around the corresponding 
main axis of the Cauchy deformation tensor C = ^ 1T ^~ 1 = RD 2 R~ l (this time both 
the state of deformation and our body rotate!). 

We see that the whole system reminds the gimbals equipped with additional mechanical device 
deforming the body. 

If our parameters D±, D 2 , Z) 3 , u are chosen in such a way that they fulfill the algebraic 
relations obtained from the equations of motion, then even when we switch off the mechanical 
device which generates the state of deformation in the material of our body, nothing will change, 
i.e., the stretchings will continue to rotate with the same constant angular velocity $ around 
the same axis in the material and the body will be rotating with the same angular velocity u 
around the same axis in the space. The deformation invariants D\, D2, -D 3 also will be constant 
during the above-described two types of rotation. 

But in our case of the affinely-rigid body subject to the Kirchhoff-Love constraints we have 
the homogeneous deformation only in the central plane of the body, whereas in the perpendicular 
direction the body performs some nonlinear oscillations. Hence, in this situation the more 
appropriate name for our special solutions is "elliptical", but we have kept the generic name 
"ellipsoidal" for the matter of convenience. 
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5 Summary 

It is interesting to note that the special solutions obtained for the polar decomposition case 
are conceptually different from those obtained for the two-polar one [6] because here the Green 
deformation tensor G = S 2 has a constant value (i.e., G = 2SS = 0) contrary to the situation 
described in [6] when the Green deformation tensor G = $ T $ = XJD 2 U~ X , as well as the Cauchy 
one C = $ _1T $ _1 = RD 2 R~ 1 , depended on time explicitly through the time dependence of U 
and R respectively, i.e., 

^ = U ($D 2 - D 2 $) U- V 0, ^ = R {uD 2 - D 2 co) RT 1 + 0, 

and performed the stationary rotations around their principal axes, whereas the deformation 
invariants A, /i, g had the constant values. 

So, if we additionally keep in mind that in [6] we obtained the stationary solutions only for 
the isotropic model J\ = J2 = J and here the general situation J\ ^ J2 is allowed, then we can 
compare the four (one in [6] and three here) studied cases according to the following scheme: 

• The only degrees of freedom we can manipulate are the rotational degrees of R- and [/-tops, 
because the deformation matrix D is constant for this type of stationary solutions. 

• To achieve the constant behaviour of the Green deformation tensor G = S 2 = (JJ = 
UD 2 U~ 1 we have to suppose that the [/-top is fixed and does not rotate at all. If U is 
constant, then the principal axes of the R- and L (= RU~ 1 )-tops (for the two-polar and 
polar decompositions respectively) rotate in the same manner, i.e., at any moment ones 
can be obtained from others with the help of applying some constant orthogonal transfor- 
mation. This situation corresponds to the above-mentioned three cases (i)-(ra) describing 
the stationary rotations of the L-top around its three principal axes. 

• If U-top is not fixed, then the Green deformation is not constant and we have to con- 
sider three branches of the stationary motion for R- and [/-tops when they rotate not 
independently but in the correlated manner, i.e., either both around their first principal 
axes or both around the second ones or both around the third ones [15]. Nevertheless, for 
our af finely-rigid body subject to the Kirchhoff-Love constraints only the third case is 
possible and exactly this situation was studied in the previous paper [6]. 

Hence, we see that in the above-described sense the results obtained in this article are essentially 
different from and simultaneously complementary to those obtained in [6]. 

Let us mention that the affine models of degrees of freedom for structured bodies have been 
studied by many authors. The thorough analysis of some stationary motions for affine bodies 
and their stability was presented in [5j EJ [TUl EH]. However, in this article we have discussed 
other problems. 
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